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摘 要 : 基于 吕梁 山北 段 管 梁山 地 区 华北 落叶 松 (Larix principis-rupprechtii ) 样 芯 资 料 , 建 立 树 轮 宽度 标准 年 表 
(STD) ,利用 Pearson 相关 分 析 方 法 得 出 3 一 4 月 平均 最 高 气温 是 研究 区 树木 径 向 生长 的 主 控 气候 因子 (R=0.509, P< 


0.05). 
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区 过 去 近 百 年 有 两 个 暧 期 和 三 个 冷 期 分 布 , 暖 基 


展 助 一 阶 线性 回归 方程 模拟 构建 近 百 年 来 管 沙 山地 区 3 一 4 月 平均 最 高 气温 的 变化 过 程 ,同时 应 用 逐一 易 
和 可 靠 性 。 对 重建 气温 序列 温度 变化 分 析 可 知 , 管 小 山地 
为 1932 


1945 年 和 1957 一 1970 年 ; 冷 期 为 1946 一 1956 年 ,1971 一 


1983 年 和 1993 一 2010 年 。 基 于 Morlet 小 波 分 析 对 重建 气温 序列 周期 分 析 可 知 ,存在 3~5 a、10~13 a、22~33 a、40~45a 
的 周期 分 布 特征 ,厄尔尼诺 -南方 涛 动 (El Nifio-Southern Oscillation ,ENSO ) 现 象 和 太阳 黑子 活动 是 其 周期 变化 驱动 


因子 。 大 尺度 空间 相关 分 析 可 知 ,重建 气温 序列 对 俄 罗 


斯 中 东部 .日 本 北部 .中国 中 东部 等 大 范围 地 区 的 温度 变化 


具有 很 好 的 空间 代表 性 。 研 究 结果 进一步 补充 和 丰富 管 小 山地 区 气候 资料 ,同时 对 森林 管理 及 农 牧 业 发 展 具有 重 


要 的 参考 意义 。 
关键 词 : 华北 落叶 松 ; 树 轮 宽 度 ， 平 均 最 高 气温 重建 ; 


气候 变化 衍生 问题 与 人 类 生存 息息相关 ,19 世 
纪 80 年 代 以 来 地 表 温度 升温 幅度 高 达 1.5 CEF, 
20 世 纪 中 期 以 后 ,大气 CO 含量 升 高 约 100ppm 之 
多 ,气候 变 暖 对 地 表 生 态 环境 影响 深远 ,对 人 类 生 
存 环境 带 来 诸多 挑战 ““。 了 解 过 去 气候 动态 进程 
对 于 系统 科学 掌握 气候 变化 规律 和 预测 未 来 气候 
变化 趋势 具有 重要 意义 。 植 被 作为 地 表 生 态 系统 
中 对 气候 响应 最 为 敏感 的 指示 物 之 一 ,可 以 精确 记 
录 和 气候 变化 的 历史 进程 。 北 半球 中 高 纬 地 区 树木 
〈 针 叶 林 树 种 ) 年 轮 具 有 分 布 范围 广 ,定年 准确 , 连 
续 性 强 ,时空 分 辩 率 高 ,样本 易 获 得 等 天 然 优势 ,成 
为 现代 许多 学 者 研究 历史 气候 变化 进程 的 主要 代 
用 指标 之 一 ”1。 

气象 站 器 测 数 据 由 于 不 同 区 域 站 点 分 布 情况 
各 异 ,气象 数据 时 间 短 ,部 分 年 限 数据 缺失 ,限制 对 
当地 长 时 间 尺 度 历史 气候 变化 的 认识 。 树 木 年 轮 
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管 小 山 


气候 学 结合 气象 站 带 测 数据 和 树 轮 宽度 数据 资料 ， 
通过 重建 历史 气候 因子 变化 情形 , 扩 宽 对 以 往 长 时 
间 尺 度 气候 变化 的 认识 。 近 年 来 树木 年 轮 气候 学 
领域 重建 历史 时 期 气温 变化 发 展 迅速 ,中 高 纬 区 域 
高 大 山脉 均 开 展 一 系列 针对 针 叶 林 树种 的 响应 和 
重建 研究 。 西 欧 阿 尔 插 斯 山 ”\ 东 亚 南部 喜马拉雅 
山西 班 牙 东 南部 内 华 达 山北 美 东部 阿 巴 拉 
契 亚 山西 亚 高 加 索 山脉 “等 地 区 , 均 开展 基于 
树木 年 轮 重建 气温 的 研究 。 在 我 国 同样 开展 大 量 
重建 历史 气温 的 研究 ,具体 区 域 分 布 于 新 疆 阿 尔 泰 
山 “\ 陕 西 秦 岭 ”、 青 海神 连 山 "河南 桐 柏 山 、 
湖北 大 别 山 全 四川 模 断 山 等 区 域 ,大 大 丰富 我 
国 基于 树木 年 轮 重建 长 时 间 尺 度 气温 研究 。 

管 水 山地 区 位 于 温带 季风 和 温带 大 陆 性 气候 ， 
半 湿 润 区 和 半 干 旱 区 分 界线 上 ,人 处 于 中 西部 地 区 气 
修 过 渡 地 带 , 当 地 树木 生长 对 气候 响应 敏感 ,适合 
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于 树木 年 轮 气候 学 研究 。 在 本 文 研究 区 域 附近 , 蔡 
秋芳 等 基于 吕梁 山中 段 油 松 样本 重建 过 去 150 a Ze 
右 5 一 7 月 平均 气温 变化 ,探究 不 同 冷 暧 期 时 段 分 区 
与 北方 大 范围 同期 气温 变化 响应 关系 ; 王 振 威 等 ” 基 
于 吕梁 山北 段 油 松 样本 重建 过 去 150 a 左右 5 一 6 月 
帕 尔 默 干旱 指数 变化 ,得 出 整个 区 域 干 旱 指 数 变化 
有 周期 和 干 湿 分 期 现象 ;李强 等 “基于 吕梁 山北 段 
油 松 样 本 重建 过 去 300 a 左右 全 年 降水 变化 ,得 出 该 
区 域 降水 年 纪 分 配 不 均匀 ,有 明显 的 周期 现象 。 在 
山西 吕梁 山区 域 树 木 年 轮 气候 学 研究 主要 集中 于 
油 松 的 研究 ,同时 针对 于 树 轮 宽度 的 响应 与 气温 重 
建 研究 比较 缺乏 。 

本 文 基于 吕梁 山北 段 管 水 山林 线 处 优势 针 叶 
林 树 种 华北 落叶 松 (Larix principis-rupprechtii ) W $È 
宽度 样 芯 资 料 ,建立 标准 年 表 , 分 析 年 表 与 各 月 不 
同时 间 段 气候 因子 的 响应 关系 ,揭示 控制 当地 树木 
生长 的 主要 气候 因子 ,重建 过 去 近 百 年 3 一 4 月 平均 
最 高 气温 变化 ,分 析 长 时 间 尺 度 下 冷暖 期 交替 、 空 
间 代 表 性 .周期 变化 等 特征 。 在 气候 变 暖 背景 下 ， 
针对 月 际 尺度 的 平均 最 高 气温 重建 研究 ,可 以 明确 
指示 历史 时 期 当地 气温 浆 值 变化 范围 ,为 应 对 未 来 
气候 变化 提供 借鉴 和 指示 。 本 文 研 究 结果 ,一定 程 
度 上 延长 需 测 气候 资料 ,丰富 我 国 中 西部 地 区 树木 
年 轮 宽度 数据 ,对 于 预测 未 来 管 洲 山地 区 气温 变化 
有 重要 作用 ,对 于 指导 当地 农林 牧 业 发 展 具有 重要 
的 参考 价值 。 
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1 材料 与 方法 


1.1 研究 区 概况 

管 沙 山地 区 隶属 于 山西 省 西部 山脉 吕梁 山北 
支 ( 图 1), 区 域内 多 年 平均 降水 量 为 390.82 mm, 多 
年 月 平均 气温 为 6.33 %C ,多 年 月 平均 相对 湿度 为 
56.18%, 月 最 低 气温 可 达 -19.06 %C, 月 最 高 气温 达 
27.25 C(I 2) ,区域 四 季 分 明 ,气温 年 差 较 大 。 夏 
季 受 暧 湿 的 东南 季风 控制 ,炎热 多 雨 ;冬季 受 干 冷 
的 西北 风 控 制 ,严寒 干燥 。 管 水 山 地 区 植被 茂密 ， 
主要 针 叶 林 树种 有 华北 落叶 松 (Larix principis-rup- 
prechtii ) , AFF (Picea meyeri) 、 油 松 (Pinus tabulaefor- 
mis) 等 ,主要 灌木 有 刺 梅 (7zorn plum) | W H (Sea 
buckthorn). RUES LL fa) eR BE MAR 
+ EAE, 
12 气象 资料 

研究 区 气象 资料 来 源 于 国家 气象 科学 数据 中 心 
(http://data.cma.cn)。 人 气象 站 点 选取 距离 采样 点 综合 
距离 最 近 的 4 个 气象 站 ,分 别 为 于 武 站 (39°00'N， 
112°18’E, 74K 1437.4 m) ,五 寨 站 (38°55'N,111°49'E， 
海拔 1401.0 m) , raf fe Hi (38°43'N, 111°35'E, rH 
1396.7 m) , AFR A (38°21'N, 111°56’E, F} 1212.2 
m)。 人 研究 中 气象 数据 具体 指标 值 均 为 4 个 气象 站 
点 数据 的 平均 值 。 该 区 域 4 个 气象 站 点 的 气候 资 
料 时 间 尺 度 为 1970 一 2020 年 ,已 有 研究 表明 该 区 域 
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图 1 管 小 山 采 样 点 及 气象 站 点 地 理 位 置 分 布 示意 图 


Fig. 1 Locations of the tree-ring sampling points and weather stations on Guancen Mountain 
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图 2 1996 一 2020 年 4 个 气象 站 多 年 月 气温 .降水 .相对 湿 
度 均 值 变化 
Fig.2 Changes in the mean monthly temperature, 
precipitation and relative humidity of four weather stations 
over the years from 1996 to 2020 


存在 气温 突变 现象 ,气温 突变 点 为 1995 一 1996 年 时 
段 内 中 。 在 确保 气候 资料 时 间 尺 度 符合 树木 年 轮 气 
候 学 研究 标准 前 提 下 ,为 保证 宽度 年 表 与 气候 因子 
相关 性 水 平 更 高 ,研究 中 选用 的 气象 数据 为 1996 一 
2020 年 的 月 降水 量 .月 平均 气温 月 最 高 气温 .月 最 
低 气温 和 月 相对 湿度 5 个 气候 指标 。 考 虑 到 树木 径 
向 生长 对 气候 因子 响应 的 “滞后 现象 ” ,树木 径 向 生长 
不 仅 受 当年 生长 季 和 气候 因 子 的 影响 ,上 一 年 生长 季 后 
期 气候 因子 同样 影响 当年 树木 的 生长 发 育 汪 7 。 基 
于 我 国 北方 针 叶 林 树种 的 生长 情况 ,选择 上 一 年 9 月 
到 当年 12 月 共计 16 个 月 的 气象 指标 进行 树 轮 宽度 
年 表 与 气候 因子 的 相关 分 析 研 究 。 
1.3 样本 采集 与 年 表 建 立 

2021 年 7 月 在 管 水 山 最 高 峰 荷 叶 坪 (2783 m) 林 
线 区 域 附近 ,选择 3 个 不 同样 点 采集 华北 落叶 松树 芯 
样本 资料 ,在 海拔 2303 m(HYPA 样 点 ) 采 集 24 棵 树 
48 根 样 芯 , 在 海拔 2477 m(HYPB 样 点 ) 和 海拔 2539 
m(HYPC 样 点 ) 分 别 采集 23 棵 树 46 根 样 芯 。 同 时 在 
管 水 山 国家 森林 公园 (HYPE 样 点 ,海拔 2128 m) 采 集 
28 棵 树 56 根 样 芯 。 采 样 中 遵循 同一 棵 树 在 胸径 和 
基 径 位 置 各 采集 一 根 样 芯 资 料 ,其 一 方便 同一 棵 树 
样本 信息 的 对 照 ,再 者 丰富 样本 量 的 同时 提升 整体 
样本 数据 的 准确 性 。 在 采样 过 程 中 ,基于 国际 树木 
年 轮 数 据 库 (International Tree- Ring Date Bank, 
ITRDB) 的 标准 ,选择 树龄 长 ,生长 健康 ,树木 长 势 好 
的 华北 落叶 松 , 利 用 内 径 为 5.15 mm 的 生长 锥 进行 样 
芯 的 采集 ,同时 确保 生长 锥 每 次 采样 均 通 过 树木 髓 
心 位 置 ,为 日 后 样 世 的 交叉 定年 ,年 表 的 制作 提供 可 
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信 度 更 高 的 数据 源 ”。 将 所 采集 的 样 世 单 根 单 简 放 
于 提前 编号 的 纸 简 内 ,防止 样 芒 资料 的 混淆 与 损 
坏 。 本 次 采样 一 共 收 集 98 棵 树 196 根 样 芯 资 料 , 采 
样 树木 胸径 在 35.62~48.90 cm。 

在 实验 室 按照 Stokes 和 Smiley 基 本 方法 对 样 世 
进行 有 序 的 拆 分 和 预 处 理工 作 2 , 样 世 安 置 于 木 覃 
中 ,依次 用 320 目 .600 目 .800 目 砂纸 进行 打磨 , 直 
至 样 忌 年 轮 年 纪 变化 清晰 可 见 ,在 双 简 显微镜 下 用 
骨架 图 法 进行 目 视 交 叉 定 年 ,确定 样 世 生长 年 限 范 
围 。 之 后 用 精度 为 0.001 mm 的 LINTAB™ 6 年 轮 宽 
度 测 量 仪 (Rinntech, Heidelberg , Germany ) 偿 年 测量 
年 轮 宽度 ,采用 COFECHA 程序 (Holmes,1983) 对 定 
年 和 测量 结果 进行 质量 检验 , 掏 除 定年 错误 或 误差 
较 大 的 样 心 ””。 最 后 使 用 ARSTAN 程序 (Cook， 
1986) 中 的 线性 函数 和 负 指 数 函 数 拟 合 树木 生长 趋 
势 , 消 除非 气候 因素 和 树木 生长 趋势 的 影响 ,通过 
双 权 重 平 均 法 进行 年 轮 曲线 的 标准 化 ” ,最终 建立 
标准 年 表 (Standard Chronology,STD ) 、 差 值 年 表 (Re- 
sidual Chronology,RES) 和 上 自 回归 年 表 (Arstan Chro- 
nology,ARS)3 种 类 型 的 年 表 。 建 立 年 表 过 程 中 , 吻 
除 掉 26 根 拟 合 树木 生长 趋势 差 的 样 芯 ,剩余 170 根 
样 必 统 一 建立 年 表 。 本 文 基于 年 表 参 数 以 及 标准 
年 表 保 留 更 多 气候 变化 低频 信和 号 信息 ,选择 标准 年 
表 (图 3) 进 行 气温 重建 分 析 *”。 

1.4 研究 方法 

利用 SPSS 软件 中 Pearson 相关 分 析 计 算 华北 落 

叶 松 树 轮 宽度 指数 年 表 与 各 月 和 特殊 月 份 组 合 气 


年 轮 宽度 指数 


1900 1920 1940 1960 1980 2000 2020 
年 份 


图 3 管 小山 华北 落叶 松树 轮 宽 度 标准 年 表 和 样本 量 


Fig.3 Standard chronology and sample sizes of the annual 


tree-ring widths of Larix principis-rupprechtii trees on 


Guancen Mountain 
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候 因 子 ( 月 平均 气温 、 月 最 低 气 温 、 月 最 高 气温 H 
降水 量 月 相对 湿度 ) 的 相关 关系 ,以 确定 影响 研究 
区 树 轮 宽度 径 向 生长 的 主要 气候 因子 。 利 用 一 阶 
线性 回归 模型 构建 气温 重建 方程 ,利用 逐一 剔除 法 
对 重建 方程 进行 交叉 检验 以 确定 方程 的 稳定 性 和 
可 靠 性 。 计 算 重建 方程 的 气温 均值 与 标准 差 ,对 
研究 区 过 去 长 时 间 气 温 变 化 进行 高 温 年 .寒冷 年 、 
偏 暖 年 . 偏 冷 年 . 冷 期 . 暖 期 划分 ”。 将 气温 大 于 
Tm+20(18.01 %C ) 的 年 份 定义 为 高 温 年 ;将 气温 介 
于 Tto(16.43 C) All Trent 20 (18.01 C) WEHE 
义 为 偏 暖 年 ;将 气温 介 于 To-a(13.28 CT) AM Tren tO 
(16.43 % ) 的 年 份 定义 为 正常 年 ;将 气温 介 于 四 
2c(11.71%) 和 Ts-c(13.28%) 的 年 份 定义 为 偏 冷 
年 ;将 气温 小 于 7,.,-20(11.71 %C ) 的 年 份 定义 为 寒 


代表 当地 华北 落叶 松 生长 情形 。 一 阶 自 相 关系 数 
H 0.683 ,表明 当年 树木 径 向 生长 受 前 一 年 气候 条 件 
的 影响 和 和 制约。 样本 总 体 解 释 量 高 达 0.936, 远 远 
超过 0.85 的 阔 值 界限 ,表明 样 芒 能 代表 该 区 域内 
华北 落叶 松树 轮 宽度 变化 的 基本 特征 。 通 过 分 析 
华北 落叶 松 样本 统计 参数 可 知 , 所 建立 的 树 轮 宽度 
标准 年 表 适 合 树 木 年 轮 气候 学 分 析 。 重 建 气温 的 
起 始 年 份 , 以 子 样本 信号 强度 SSS(Sub-sample Sig- 
nal Strength ) 来 确定 ,为 尽 可 能 的 延长 重建 气温 时 
间 长 度 , 选 取 SSS BEX 0.750% ,结合 年 表 参 数值 ， 
得 出 可 信 年 表 时 间 尺 度 为 1923 一 2020 年 ( 表 1)。 
2.2 树木 径 向 生长 与 气候 因子 的 相关 关系 

考虑 到 气候 因子 对 树木 径 向 生长 的 滞后 效应 ， 
选择 上 一 年 9 月 至 当年 12 月 的 月 平均 气温 .月 最 低 


冷 年 ”“”"。 为 进一步 分 析 管 小 山地 区 3 一 4 月 平均 最 
高 气温 的 年 代 际 变化 特征 ,对 重建 气温 序列 进行 11 
a 滑动 平均 ,按照 气温 重建 序列 连续 11 a 大 于 或 等 
F Tacan (14.86 C) 划 分 为 暧 期 ,连续 11 a NF Tran RI 
分 为 冷 HH), ALA Matlab 软件 对 重建 气温 序列 进 
行 Morlet 小 波 分 析 , 以 探究 气温 变化 是 否 存 在 周期 现 
象 ”“。 利 用 荷兰 星 家 气象 研究 所 数据 共享 网 站 
(http://www.climexp.knmi.nl) ,对 气温 重建 序列 进行 
空间 代表 性 分 析 。 


2 结果 与 分 析 
2.1 年 表 特 征 分 析 


气温 月 最 高 气温 .月 降水 量 月 相对 湿度 气候 因子 
与 树 轮 宽度 标准 年 表 做 Pearson 相关 分 析 ( 图 4)。 
结果 表明 , 树 轮 宽 度 指数 年 表 与 绝 大 部 分 月 份 降 水 
量 相 关 性 不 显著 ,只 与 当年 3 月 降水 量 呈 显著 负 相 
关 关 系 (R=-0.364,P<0.05)。 树 轮 宽度 指数 年 表 与 
当年 1 月 (R=0.429, P<0.05) 、3 月 (R=0.537, P< 
0.05) 、4 H (R=0.447, P<0.05) 、5 月 (R=0.407, P< 
0.05) .10 月 (R=0.446,P<0.05) 的 月 最 高 气温 呈 显 著 
正 相 关 关 系 , 与 3 月 平均 气温 (R=0.433,P<0.05) 和 6 
月 最 低 气 温 (R=0.414,P<0.05) 明 显著 正 相 关 关 系 。 
Pearson 相关 性 分 析 可 知 , 管 水 山地 区 华北 落叶 松 径 
向 生长 主要 受 月 最 高 气温 限制 ,通过 组 合 月 份 相 关 
性 分 析 可 知 , 树 轮 宽度 指数 年 表 与 3 一 4 月 平均 最 高 


由 采样 点 华北 落叶 松树 轮 宽 度 标 准 年 表 信息 
可 知 ,平均 敏感 度 为 0.213 , 信 品 比 高 达 30.299 ,说 明 
采样 点 标准 年 表 保留 有 丰富 的 气候 信息 。 标 准 差 
为 0.302 ,序列 间 .树木 内 、 树 木 间 有 较 高 的 平均 相关 
系数 ,表明 采样 点 年 表 有 较 好 的 区 域 一 致 性 ,可 以 


表 1 


气温 呈 显 著 正 相关 关系 (R=0.509,P<0.05) ,3 一 4 月 

平均 最 高 气温 为 限制 当地 优势 种 华北 落叶 松 径 向 

生长 的 最 主要 气候 因子 。 

2.3 过 去 近 百 年 3 一 4 月 平均 最 高 气温 重建 与 检验 
基于 Pearson 相关 分 析 结 果 ,将 管 小 山地 区 华北 


管 小山 华北 落叶 松 采样 点 概况 及 树 轮 样本 统计 参数 


Tab.1 Sampling points and statistical parameters of tree-ring samples of Larix principis-rupprechtii 


trees on Guancen Mountain 


特征 参数 统计 值 特征 参数 统计 值 

采样 点 海拔 2128~2539 m 序列 间 的 平均 相关 系数 (R) 0.490 

样本 量 ( 树 / 心 ) 98/196 树木 内 的 平均 相关 系数 (R) 0.514 

序列 长 度 1899 一 2020 年 树木 间 的 平均 相关 系数 (R;) 0.482 

公共 区 间 1955 一 2020 年 信 品 比 (SNR) 30.299 
平均 敏感 度 (MS) 0.213 样本 总 体 解释 量 (EPS) 0.936 

标准 差 (SD) 0.302 一 阶 自 相关 系数 (4C') 0.683 

子 样本 信号 强度 SSS > 0.75 起 始 年 ( 样 世 数 ) ”1923 年 (6) 
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O 月 平均 气温 Œ 月 最 低 气 温 GS 月 最 高 气温 
BAe Bw 月 相对 湿度 


a=0.05 


相关 系数 


-0.3 a=0.05 


注 :p 代 表 前 一 年 ;ic 代表 当 年 。 
图 4 管 水 山 华 北 落叶 松树 轮 宽 度 标 准 年 表 与 气候 因子 的 
相关 分 析 


Fig. 4 Correlation analysis of standard chronology of tree- 


ring width and climate factors of Larix principis-rupprechtii 


trees on Guancen Mountain 


落叶 松树 轮 宽度 标准 年 表 值 作为 自 变 量 ,3 一 4 月 平 
均 最 高 气温 作为 因 变 量 , 采 用 线性 回归 方法 构建 管 
深山 地 区 气温 重建 方程 : 
T a34 = 7.4074GCsw + 7.5207 

式 中 :7 为 3 一 4 月 平均 最 高 气温 ,CCsm 为 管 小 山 
华北 落叶 松树 轮 宽度 标准 年 表 值 。 

该 重建 方程 的 相关 系数 为 0.51, 可 以 解释 该 区 
域 3 一 4 月 平均 最 高 气温 序列 方差 ( 尼 ) 的 25.89% , 调 
整 自由 度 后 的 方差 解释 量 ( 尼 ) 为 22.67% ,下 检验 值 
为 8.034,P<0.01。 

重建 方程 的 稳定 性 和 可 靠 性 直接 影响 到 重建 
序列 的 质量 ,采用 逐一 剔除 法 对 气温 重建 方程 进行 
交叉 检验 ,验证 重建 方程 的 稳定 性 和 可 靠 性 ,各 项 
检验 统计 量 数值 见 表 2。 符 号 检验 (0S7) 和 威 尔 科 交 
森 符号 秩 检 验 (Wilcoxon S7T) 均 达到 了 显著 水 平 (P< 
0.01) ,说 明 3 一 4 月 平均 最 高 气温 重建 序列 与 实测 
序列 在 高 低频 变化 上 有 较 好 的 一 致 性 ;误差 缩减 值 
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(RE) 是 精确 检验 估计 气候 重建 值 可 靠 性 的 统计 量 ， 
取 值 范围 为 L[- ,1] ,数值 小 于 或 等 于 0 表示 重建 结 
果 不 可 信 , 大 于 0 表示 重建 结果 具有 一 定 可 靠 性 , 越 
接近 1 表示 重建 结果 越 可 靠 ” ,本 文 RE 为 0.95 ,说明 
气候 重建 值 通 过 检验 ;乘积 平均 数 (1) 可 以 双重 考虑 
到 估计 预报 量 和 实际 预报 量 的 符号 及 数量 值 , 本 文 ; 
值 为 3.00, 通 过 显著 性 检验 (P<0.01 )。 基 于 逐一 史 
除法 对 重建 方程 的 交叉 检验 分 析 ,重建 管 水 山地 区 
3 一 4 月 平均 最 高 气温 方程 通过 了 稳定 性 和 可 靠 性 
检验 。 从 图 5 可知, 重建 气温 与 实测 气温 值 有 很 好 
的 同步 性 。 
2.4 重建 气温 序列 波动 变化 特征 

由 网 6 可 知 , 过 去 近 百 年 来 , 管 水 山地 区 3 一 4 
月 平均 最 高 气温 动态 波动 明显 ,20 世纪 末期 以 来 ， 
气温 整体 旦 波动 上 升 趋势 。 重 建 的 3 一 4 月 平均 最 
高 气温 波动 范围 在 10.94 % (1996 年 ) ~18.29 C 
(2015 年 ) ,整体 气温 相差 为 7.35 Co TEES BE 
段 ,3 一 4 月 平均 最 高 气温 的 平均 值 (7,) 为 
14.86 ,标准 差 (o) 为 1.57 Co 

由 图 6 可 知 ,在 过 去 近 百 年 里 , 管 沙 山地 区 高 温 
年 出 现 1 次 ,为 2015 年 , 占 总 年 份 的 1.02% ;寒冷 年 
出 现 4 次 ,分 别 为 1923 年 .1996 年 .1997 年 .1998 年 ， 
占 总 年 份 的 4.08% 5 偏 暖 年 与 偏 冷 年 均 出 现 13 次 ， 
分 别 占 总 年 份 的 13.27%; 正 常年 为 67 次 , 占 总 年 份 
的 68.36%。 整 个 长 时 间 尺 度 重建 序列 ,正常 年 居 
多 ,表明 管 小 山地 区 3 一 4 月 平均 最 高 气温 波动 变化 
稳定 。 

在 过 去 近 百 年 来 , 管 沙 山地 区 经 历 了 两 个 暖 期 
和 三 个 冷 期 ( 表 3)。 两 个 暖 期 持续 时 间 均 为 14 a, 后 
一 个 暖 期 1957 一 1970 年 气温 整体 明显 高 于 前 一 个 
BHH 1932—1945 年 ,后 一 个 暖 期 期 间 气 温 更 高 ,更 
炎热 。 冷 期 持续 时 间 最 长 (1993 一 2010 年 ) 为 18 a, 
同时 也 是 冷 期 中 气温 最 寒冷 的 阶段 ,气温 距 平 值 
为 -0.643 %C; 冷 期 持续 时 间 最 短 (1946 一 1956 年 ) 为 
11 a, 是 冷 期 中 相对 温暖 的 阶段 ,气温 距 平 值 为 
-0.169 ;另外 一 个 冷 期 (1971 一 1983 年 ) 持 续 时 间 


表 2 逐一 剔除 法 检验 统计 量 


Tab. 2 Statistical characteristics of leave-one-out test 


时 段 r R Ra 


F ST Wilcoxon ST RE t 


1996—2020 4E 0.51 25.89% 22.67% 


RE 为 误差 缩减 值 ;i 为 乘积 平均 数 ;** 代 表 99% 的 置信 水 平 。 


8.04” 16+/9-~ 16+/9-~ 0.95 3.00” 


注 :为 相关 系数 ; 尼 为 方差 解释 量 ; 忆 为 调整 自由 度 后 的 方差 解释 量 ;为 方差 检验 值 ;$S7 为 符号 检验 ;Wilecoxon ST 为 威 尔 科 克 森 符 号 秩 检 验 ; 
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| pel 为 13 a, 气 温 距 平 值 为 -0.437 %。 对 管 小 山地 区 3 一 
4 月 平均 最 高 气温 年 代 际 变化 特征 分 析 可 知 ,不 同 

24 冷暖 期 旦 波动 变化 特征 ,冷暖 期 交替 分 布 。 

m 25 重建 气温 序列 周期 变化 特征 

S 使 用 Morlet 连续 小 波 变 换 对 气温 重建 序列 进行 
mg 16 周期 分 析 ,为 减弱 连续 小 波 变 换 对 始 ,未 区 域 的 边界 
本 效应 ,选择 对 称 延 伸 法 来 减弱 边界 效应 的 影响 “2 。 
由 网 7 可 知 , 管 沙 山地 区 近 百 年 3 一 4 月 平均 最 高 气 
温 重建 序列 有 明显 的 周期 变化 特征 。 在 0~30 a 主 
gle AT BEA 3 ~5 a 准 周期 变化 存在 于 整个 序列 中 ; 
SEP FH HF sh gd GP Pg 1970 年 之 后 ,存在 明显 的 10~13 a 准 周期 变化 特 
年 份 征 。 在 30~60 a 主 周期 阶段 内 ,22~33 a 和 40~45 a 准 


周期 变化 存在 于 整个 序列 中 。 
2.6 重建 气温 序列 空间 代表 性 分 析 


图 5 管 水 山地 区 3 一 4 月 平均 最 高 气温 重建 值 与 实测 值 对 比 


Fig. 5 Comparison between the reconstructed value of the 


Moth EER BS Ate yy Ve. fey Se er 
mean maximum temperature from March to April and the 为 进一步 探究 管 沙 山地 区 3 一 4 月 平均 最 高 气 
measured value in Guancen Mountain area 温 重 建 序列 能 和 否 代 表 大 范围 的 温度 变化 ,本 研究 利 

| —— 重建 气温 一 一 11 a 滑动 平均 什 


温度 /C 


1930 1940 1950 1960 1970 1980 1990 2000 2010 2020 
年 份 


图 6 管 水 山地 区 3 一 4 月 平均 最 高 气温 重建 序列 及 11 a 滑动 平均 值 


Fig.6 The reconstruction sequence of the mean maximum temperature from March to April and the 11 year moving 


average value in Guancen Mountain area 


表 3 重建 气温 序列 冷暖 期 划分 


Tab. 3 Reconstructed temperature Series division of cold and warm periods 


暧 期 年 数 /a PROF A/C 冷 期 年 数 /a 距 平 值 /*C 
1932 一 1945 年 14 +0.249 1946 一 1956 年 11 -0.169 
1957—1970 4 14 +0.589 1971—1983 4E 13 -0.437 


1993—20104F 18 -0.643 
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图 7 管 深 山地 区 3 一 4 月 平均 最 高 气温 重建 序列 
的 小 波 分 析 
Fig. 7 Wavelet analysis of reconstructed temperature 
sequence of the mean maximum temperature from March to 


April in Guancen Mountain area 


用 咒 测 和 重建 气温 资料 分 别 与 同时 间 段 的 气温 格 
点 数据 (CRU TS 4.05 ,0.5°x0.5° ) 做 空间 相关 性 分 析 
(图 8)。 结 果 表 明 ,器 测 与 重建 3 一 4 月 平均 最 高 气 
温 与 格 点 数据 空间 相关 模 态 基本 一 致 ,两 者 相关 性 
最 好 的 区 域 (R>0.4) 均 为 山西 北部 地 区 。 同 时 重建 
的 管 小 山地 区 3 一 4 月 平均 最 高 气温 序列 也 能 够 较 
好 的 代表 俄罗斯 中 东部 .日 本 北部 .中国 中 东部 等 
大 范围 地 区 的 温度 变化 。 通 过 对 比 器 测 和 重建 数 
据 与 气温 格 点 数据 的 空间 相关 性 分 析 可 知 ,器 测 数 
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据 的 空间 代表 性 范围 更 大 ,相关 程度 更 高 ,因此 在 
未 来 该 区 域 基于 树 轮 宽度 重建 古 气 候 的 研究 中 , 需 
要 寻找 年 龄 跨度 更 大 ,对 气候 更 敏感 的 树 轮 宽 度 样 
本 资料 。 


3 讨论 


31 管 水 山 华北 落叶 松 径 向 生长 对 气候 因子 的 
响应 

华北 落叶 松树 轮 宽 度 指 数 年 表 与 3 月 降水 量 呈 
显著 负 相 关 关 系 , 管 小 山地 区 属于 温带 季风 气候 和 
温带 大 陆 性 气候 过 渡 地 带 ,3 月 雨季 未 来 临 ,本 身 区 
域 降水 较 少 ,但 3 月 气温 回暖 速度 较 快 ,河流 解冻 以 
及 当地 地 下 水 资源 补充 树木 生长 所 需 。 早 春 树木 
生长 发 育 需 水 量 大 ,但 存在 阔 值 范围 ,过 量 的 水 分 
条 件 抑 制 树木 光合 作用 和 呼吸 作用 的 进行 ,阻碍 树 
木 体 内 能 量 物质 的 传输 ”。 树 轮 宽度 指数 年 表 与 1 
月 最 高 气温 晨 显 著 正 相关 关系 ,在 冬季 ,呼吸 作用 
是 树木 最 主要 的 生理 活动 ,温暖 的 冬季 减少 树木 体 
内 储存 的 碳水 化 合 物 的 流失 ,更 有 利于 下 一 年 树木 
的 生长 发 育 ;相反 严寒 的 冬季 ,温度 过 低 , 对 树木 体 
内 的 植物 组 织 和 细胞 造成 永久 性 损伤 ,抑制 树木 生 
长 发 育 ”。 树 轮 宽度 指 数 年 表 与 春季 月 份 最 高 气 
温 呈 显著 正 相 关 关系 ,春季 是 树木 生理 活动 和 径 问 
生长 的 旺盛 时 期 ,温度 的 升 高 ,促进 华北 落叶 松 形 
成 层 活动 的 提早 开始 ,同时 加 速 积 雪 融 化 ,利于 当 
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图 8 管 小 山地 区 1996 一 2020 年 3 一 4 月 平均 最 高 气温 器 测 数 据 和 重建 数据 与 CRU 格 点 数据 的 空间 相关 分 析 
Fig. 8 Spatial correlation analysis of the measured and reconstructed data of the mean maximum temperature from March to April 
1996 to 2020 and the CRU grid data in Guancen Mountain area 
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年 树木 生长 发 育 , 促 进 宽 轮 的 形成 3。 同时 张 
文 涛 等 "在 对 芦 芽 山 阳 坡 华北 落叶 松 和 白 杆 的 研 
究 中 , 均 发 现 春季 气温 促进 树木 生长 发 育 的 结论 ， 
进一步 验证 本 文 研究 结果 的 合理 性 。 树 轮 宽 度 指 
数 年 表 与 6 月 最 低 气 温 呈 显著 正 相 关 关 系 ,6 月 为 春 
夏 交 接 的 季节 ,气温 波动 幅度 较 大 ,保持 树木 生长 
底 温 的 需要 ,对 于 维持 轮 宽 的 发 育 至 关 重 要 。 树 轮 
宽度 指数 年 表 与 10 月 最 高 气温 呈 显 著 正 相关 关系 ， 
10 月 处 于 华北 落叶 松 生长 季 末 期 ,树木 并 未 休 眼 ， 
较 高 的 温度 延长 树木 生长 季 ,利于 光合 作物 的 积 
累 ,对 晚 材 发 育 起 促进 作用 3H。 本 文 基于 Pearson 相 
关 分 析 结 果 和 树木 生理 学 原理 ,选择 对 管 小 山地 区 
华北 落叶 松 生长 影响 最 为 显著 的 3 一 4 月 平均 最 高 
气温 气候 因子 进行 近 百 年 气温 重建 研究 。 
3.2 重建 气温 序列 对 比 与 分 析 

为 了 进一步 验证 管 小 山地 区 近 百 年 来 3 一 4 月 
平均 最 高 气温 重建 结果 的 准确 性 和 可 靠 性 ,将 本 文 
重建 气温 序列 与 陕西 镇 安县 3 一 4 月 平均 最 高 气温 重 
EPI! 秦岭 太白 山 3 一 6 月 平均 气温 重建 序列 
做 对 比分 析 (图 9)。 陕 西 镇 安县 重建 气温 序列 
1932—1937 4 5j 1941—1945 年 两 个 暖 期 和 1953 一 
1956 年 冷 期 与 本 文 重建 结果 相 一 致 ;秦岭 太白 山 重 
建 气温 序列 1977 一 1983 年 和 1993 一 1996 年 两 个 冷 
期 与 本 文 重建 结果 相 一 致 。 本 文 的 寒冷 年 (1996 一 


温度 /'C 
加 
a 


(a) 本 文 管 小 山地 区 3 一 4 月 平均 最 高 气温 A 


1998 年 ) 分 布 时 间 段 与 陕西 镇 安县 和 太白 山 研究 结 
果 基 本 相同 。 本 文 唯一 的 高 温 年 为 2015 年 ,该 年 厄 
尔 尼 诺 现象 是 有 观测 记录 以 来 最 强 的 3 次 厄尔尼诺 
现象 之 一 ,同时 也 是 维持 时 间 最 长 的 厄尔尼诺 过 
程 ,造成 北方 旱 区 伏 秋 旱 严 重 , 农 业 和 畜牧 业 受 到 
不 同 程度 的 影响 汪汪 ,进一步 验证 本 文 重建 结果 的 
合理 性 。 本 文 气温 重建 结果 与 周边 气温 重建 情形 
基本 一 致 ,说 明 在 大 范围 空间 尺度 上 ,重建 管 小 山 
过 去 近 百 年 3 一 4 月 平均 最 高 气温 变化 有 较 好 的 代 
表 性 。 
3.3 重建 周期 与 内 在 驱动 力 

本 文 气温 重建 序列 存在 明显 的 3~5 a、10~13 a、 
22~33 a、40~45 a 的 准 周期 变化 特征 。 厄 和 尔 尼 诺 - 南 
方 涛 动 (El Nifio-Southern Oscillation ,ENSO ) 现 象 的 
周期 为 2~7 a, 与 本 文 3~5a 准 周期 变化 相 一 致 ,说 明 
ENSO 现象 对 管 水 山 树木 径 向 生长 影响 深远 ,该 地 
区 气温 变化 与 大 尺度 海 气 间 耦 合 振荡 有 关 ,ENSO 是 
影响 该 区 域 气温 变化 的 重要 因素 之 一 “2 。 太 阳 黑 
子 活动 的 周期 在 11 a 左右 ,与 本 文 10~13a 和 11a 售 
数 的 周期 相 一 致 ,进一步 说 明太 阳 黑 子 活 动 对 管 小 
山区 域 温 度 变 化 有 重要 影响 所 。 在 我 国 河南 尧 山 汪 、 
新 疆 阿 尔 泰山 中 湖北 神农 架 地 区 秦岭 太白 山 呈 、 
青海 祁连山 "J 川 汗 横 断 山 中 和 吉林 长 白山 中 气温 
重建 过 程 中 , 均 发 现 ENSO 现象 或 者 太阳 黑子 活动 
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图 9 管 小山 重建 气温 与 其 他 区 域 重 建 气温 序列 对 比 


Fig. 9 Comparison of the reconstructed temperature in Guancen Mountain and other regions 
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对 温度 变化 影响 深远 。 通 过 其 他 区 域 气 温 重 建 周 
期 的 认识 ,进一步 说 明 本 文 重 建 气 温 中 ,周期 存在 
的 合理 性 ,这 些 周 期 现象 的 存在 ,暗示 管 小 山地 区 
与 全 球 气候 变化 的 内 在 关联 ,以 及 受到 大 范围 气候 
变化 的 影响 显著 。 


4 结论 


利用 采样 自 吕 深山 北 受 管 水 山区 的 华北 落叶 
松树 轮 宽度 样 芯 资料 ,建立 该 区 域 树 轮 宽度 指数 标 
准 年 表 ， nate ara eared ae 
系 , 进 而 重建 该 区 域 近 百年 来 3 一 4 月 平均 最 高 气温 

的 变化 ,经 过 对 重建 气温 序列 检验 Re Jal 

期 分 析 和 空间 代表 性 分 析 ,得 到 以 下 结 

(1) 3 一 4 月 平均 最 高 气 sarees ieee 
北 落叶 松 径 向 生长 最 主要 的 气候 因子 。 

(2) 重建 结果 可 知 , 管 水 山地 区 过 去 近 百 年 有 
两 个 暖 期 和 三 个 冷 期 分 布 , 暖 期 为 1932 一 1945 年 和 
1957 一 1970 年 ; 冷 期 为 1946 一 1956 年 ,1971 一 1983 
年 和 1993 一 2010 年 。 

(3) 重建 气温 序列 存在 3~5 a, 10~13 a, 22~33 
a, 40~45 a 的 周期 分 布 特征 ,厄尔尼诺 -南方 涛 动 
(ENSO) 现 象 和 太阳 黑子 活动 是 其 周期 变化 驱动 
因子 。 

(4) 重建 气温 序列 对 俄罗斯 中 东部 .日 本 北部 、 
中 国 中 东部 等 大 范围 地 区 的 温度 变化 具有 很 好 的 
空间 代表 性 
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Variation of mean maximum temperature from March to April since 1923 in the 


northern section of Lvliang Mountain recorded by tree rings 


WANG Jiachuan'”, LI Shuheng™, GUO Yili, HAN Yijie”, MAO Zhonglei” 
(1. College of Urban and Environmental Sciences, Northwest University, Xi’ an 710127, Shaanxi, China; 2. Shaanxi 
Key Laboratory of Earth Surface System and Environmental Carrying Capacity, Xi an 710127, Shaanxi, China) 


Abstract: With the established standard chronology of the tree-ring width of Larix principis-rupprechtii in the 
Guancen Mountain in the northern section of Lvliang Mountain as basis, Pearson correlation analysis concluded 
that the average maximum temperature from March to April was the main controlling climate factor for the radial 
growth of trees in the study area (R = 0.509, P < 0.05). First-order linear regression equation was used to simulate 
and construct the change process of the average maximum temperature in the Guancen Mountains from March to 
April in the past nearly a century. The reconstructed equation was then cross-checked by the elimination method 
of leave-one-out test to verify its stability and reliability. Analysis of the temperature changes of the reconstructed 
temperature series showed that two warm periods and three cold periods in the past nearly a century in the 
Guancen Mountain area; the warm periods were 1932-1945 and 1957-1970, and the cold periods were 1946- 
1956, 1971-1983, and 1993-2010. Periodic analysis of the reconstructed temperature series based on Morlet 
wavelet analysis revealed periodic distribution characteristics of 3~5 a, 10~13 a, 22~33 a, and 40~45 a. El Niño- 
Southern Oscillation and sunspot activity were the drivers of this periodic variation. Large- scale spatial 
correlation analysis indicated that the reconstructed temperature series had a good spatial representation of 
temperature changes in large-scale regions such as central and eastern Russia, northern Japan, and central and 
eastern China. These results have important reference value for forest management and agriculture and animal 
husbandry development in Guancen Mountain area. 

Keywords: Larix principis-rupprechtii; tree-ring width; mean maximum temperature reconstruction; Guancen 


Mountain 


